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Abstract. Dynamic behavior of a site diluted Ising ferromagnet in the presence of 

periodically oscillating magnetic field has been analyzed by means of the effective 

field theory (EFT). Dynamic equation of motion have been solved for a honeycomb 

tt \ lattice (z — 3) with the help of a Glauber type stochastic process. The global phase 

diagrams and the variation of the corresponding dynamic order parameter as a function 

of the Hamiltonian parameters and temperature has been investigated in detail and it 

has been shown that the system exhibits reentrant phenomena, as well as a dynamic 

tricritical point which disappears for sufficiently weak dilution. 
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1. Introduction 

Ising model under a time dependent external magnetic field (kinetic Ising model) has 
attracted much interest from both theoretical and experimental points of view recently. 
Due to the competing time scales of the system and periodic external magnetic field, 
the system cannot respond to the external magnetic field instantaneously which causes 
interesting behaviors. At high temperatures for the high amplitudes of the periodic 
magnetic field, the system follows the external field with a delay while this is not 
the case for low temperatures and small magnetic field amplitudes. This spontaneous 
symmetry breaking indicates the presence of a dynamical phase transition (DPT)[T] 
which shows itself in the dynamical order parameter (DOP) which is defined as the time 
average of the magnetization over a full period of the oscillating field. DPT is firstly 
observed theoretically within the mean field approximation (MFA) [2J. Since then, much 
attention has been devoted to investigate the dynamic nature of the phase transition 
in experimental studies. As an example, a DPT occurring in high frequency region is 
represented by Jiang et al. [3j using the surface magneto-optical Kerr-effect technique 
for epitaxially grown ultrathin Co films on a Cu(OOl) surface. For a [Co(4A°)/Pt(7A°)] 
multilayer system with strong perpendicular anisotropy, an example of DPT is shown 
by Robb et al. [I]. In this study, they found that the experimental nonequilibrium 
phase diagrams are found to strongly resemble the dynamic behavior calculated from 
simulations of a kinetic Ising model. Thus, a strong evidence of consistency between 
theoretical and experimental studies is shown. Polyethylene naphtalate nanocomposites 
are investigated by Kanuga et al. [5]. After detailed analysis, they shown that the 
material undergoes three critical structural transitions. 

DPT and hysteresis behaviors of the S-l/2 Ising model with periodic external 
magnetic field has been studied by means of MFA [6j [7J [8j El [10] and Monte Carlo 
(MC) simulations [HJ [121 Q31 CHI EEl EEl E], and there are also some works for higher 
spins [HI [19] or mixed spin Ising systems [201 [21] , as well as the other variants of the 
model such as the transverse kinetic Ising model [22J based on the MFA. 

It is well known that MFA neglects the thermal fluctuations via neglecting the 
self spin correlations whereas the effective field theory (EFT) which is superior to 
conventional MFA is based on the differential operator technique [231 EI] and it takes 
into account the self spin correlations. Thus, it is expected from EFT to obtain more 
reasonable results than MFA for these systems, as in the case of static Ising model. EFT 
has been applied to the kinetic Ising model in the past few years. For example, thermal 
and magnetic properties of kinetic Ising model [25], EH ETJ EH], and also the transverse 
kinetic Ising model [291 30j have been studied by EFT very recently. Based on the results 
obtained in these works, we can conclude that EFT gives more reasonable results for the 
investigation of the DPT behaviors of the kinetic Ising model, in comparison to MFA. 

On the other hand, as far as we know, the thermal and magnetic properties of 
the quenched disordered systems such as site (or bond) diluted ferromagnets driven 
by a periodically oscillating magnetic field have not yet been investigated within the 



framework of EFT or MFA . Furthermore, DPT properties of such disordered materials 
have not yet also been considered by means of MC simulations. These types of disorder 
effects constitute an important role in material science, since the quenched disorder 
effects may induce some important macroscopic effects on the material. Therefore, we 
believe that the investigation of quenched randomness effects on the DPT properties 
of Ising model still needs particular attention and we need a detailed consideration 
of competition between the disorder conditions and the oscillation of magnetic field. 
Hence, in this work we intend to study the DPT properties of a kinetic Ising system by 
introducing the quenched site dilution effects. For this purpose, we organized the paper 
as follows: In Sec. [2] we briefly present the formulations. The results and discussions 
are summarized in Sec. |3l and finally Sec. H] contains our conclusions. 

2. Formulation 

We consider an Ising ferromagnet defined on a lattice which has a coordination number 
z with a time dependent external magnetic field. The Hamiltonian describing our model 

is 

H = -J 5Z c i°j s i s j - H(t)^2dSi, (1) 

where J > is spin-spin exchange interaction, q is a site occupation variable and Sj is 
the spin variable. Site occupation variable can take the values c$ = which means that 
the site i is empty or q = 1 if the site i is occupied by a magnetic atom and the spin 
variable can take values Si = ±1. The first summation in Eq. (OQ) is over the nearest 
neighbor site pairs and the second one is over all lattice sites. The time dependent 
external magnetic field is given by 

H{t) = H cos(ut), (2) 

where t is the time and H is the amplitude of the oscillating magnetic field with 
a frequency u>. The dynamical evolution of the system may be given by Glauber 
dynamics[31j based on a master equation 

Tj{{c iSi )) r = -{{c iSi )) r + ((a tanh \Pa (E t + H(t))])) r , (3) 

where 1/r is the transition per unit time in a Glauber type stochastic process, 
(3 = l/ksT and ks represents the Boltzmann constant, T is the temperature and Ei is 
the local field acting on the site i and it is given by 

z 

Ei = jJ2 c 6S6, (4) 

5=1 

where z is the coordination number of the lattice. The inner average brackets in Eq. ([3]) 
stands for the thermal average and the outer one (which has a subscript r) represents 
the random configurational average which is necessary for including the site dilution 

effects. 



In order to handle the second term on the right hand side of the Eq. ([3]) one can use 
the differential operator technique [23, 24J. By using the differential operator technique, 
Eq. (EJ) gets the form 

r—((ciSi)) r = -((ciSi)) r + ((ciexp(ciEiV))) r f(x), (5) 

where V = d/dx is the differential operator and the function f(x) is given by 

f(x)=tanh[/3(x + H(t))}. (6) 

The effect of the differential operator on a function f(x) is given by 

exp(aV)/(x) = /(x + o), (7) 

with any constant a. By taking into account the two possible values of q as q = 0, 1 
we can write the exponential term as 

exp (aci) = Ci exp (a) + f — Cj, (8) 

where a is any constant. By using Eq. ([8]) in Eq. (jSJ) we can obtain 

T-T,{{CiSi)) r = -((CiSi)) r + ((Cj exp(EjV))) r /(x), (9) 

where for site occupations c] = Ci was used. By using Eq. (TJJ in Eq. fl9]) we get 

r-7-({ciSi)) r = -((c i s i )) r + ( (a Y[ [c s exp(Js s V) + 1 - c 5 ] \\ f(x). (10) 

In order to get a polynomial form of the second term on the right hand side of Eq. f JT0|) , 
we write the exponential term in terms of the hyperbolic trigonometric functions 

T -7;(( c i s i))r = -((CiSi))r 

II * \\ (11) 

+ ((ciY[[c s cosh( JV) + c s ss sinh( JV) + l-c s }j) f(x). 

When the product in Eq. ffTTl) is expanded, the multi site spin correlations appear. For 
simplicity let us handle these correlations with an improved decoupling approximation 
(DA) [32] as 

((ci . . . CjC k s k . . . cisi)) r = ((c i )) r ...((c j )) r ((c k s k )) r ...((cisi)) r , (12) 

with 

(ci) r = c, ((ciSi)) r = m. (13) 

Detailed discussion of DA in site dilution problem for the static case can be found in[32j. 
By using Eqs. (fT2|) and (TTBl in Eq. [TT] we get 

(XTfl 

t— = -m + c(([ccosh(JV)+msinh(JV) + I-c] 2 )) r /(x), (14) 

By using the binomial expansion and writing the hyperbolic trigonometric functions in 
terms of the exponential functions we get the compact form of Eq. f !14p as 

dm 1 / ^ \ 

-— = - -m + > A r m r , (15) 

dt r\ ^ J 



where 

Ar = i2(j( P ) c p - r+1 (l - c) z - p cosh p - r (JV) sinh r (JV)/(x). (16) 

Eq. f|T5l) can be solved by various numerical methods. In this work we prefer to use 
the fourth order Runge-Kutta method (RK4) by regarding Eq. (TT5|) as an initial value 
problem. We can mention that the differential equation derived in Eq. (fT5l) extends up 
to the term m z . Each term in Eq. (TT5]) makes contribution to the solution, because the 
value of m, which is calculated at each time step, is related to the previous m value 
however the situation is different from the behavior of the equilibrium systems at which 
high ordered terms can be neglected in the neighborhood of phase transition point. 

The system has three Hamiltonian parameters, namely magnetic field frequency 
(cu), magnetic field amplitude (H ) and site concentration (c). For selected values of 
these Hamiltonian parameters and temperature, RK4 will give convergency behavior 
after some iterations i.e. the solutions have property m(t) = m(t + 2tt/u) for arbitrary 
initial value for the magnetization (m(t = 0)). Thus, after obtaining this convergent 
region after some transient steps (which depends on Hamiltonian parameters and the 
temperature) the DOP can be calculated from 

Q = y I m{t)dt. (17) 

The dynamic nature of the system for a given temperature and Hamiltonian 
parameters can be determined from the DOP. There are three possible states for the 
system, namely ferromagnetic (F), paramagnetic (P) and mixed (F+P). The solution 
m(t) in the convergent region satisfies 

m(t) = -m(t + n/u), (18) 

in the P phase which is called the symmetric solution. The solution corresponding 
to P phase follows the external magnetic field and oscillates around zero value which 
means that the DOP Q is zero. In the F phase, the solution does not satisfy Eq. ( fi~8l 
and this solution is called as non-symmetric solution which oscillates around a non-zero 
magnetization value, and does not follow the external magnetic field i.e. the value of Q 
is different from zero. In these two cases, the observed behavior of the magnetization is 
regardless of the choice of the initial value of magnetization m(0) whereas the last phase 
has magnetization solutions symmetric or non-symmetric depending on the choice of 
the initial value of magnetization corresponding to the coexistence region where F and 
P phases overlap. 

3. Results and Discussion 

In order to determine the effect of the Hamiltonian parameters on the critical 
temperature of the system we present the phase diagrams for a honeycomb lattice (z = 3) 
in (kBT c /J — Ho/ J) and (IcbT c /J — c) planes with frequencies to = 0.1,0.5, 1.0,5.0 in 
Figs. (01 and (@J, respectively. 
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At first sight we can see from Fig. (J3J) that as c value decreases then the F phase 
region gets narrower in the (kBT c /J — Ho/ J) plane. This is an expected result, since as 
the lattice sites are diluted more and more then the energy contribution which comes 
from the spin-spin interactions gets smaller. Hence, the system can undergo a DPT 
at smaller critical temperatures, due to the energy originating from the temperature 
and (or) magnetic field which overcomes the spin-spin interactions. Therefore, the 
ferromagnetic regions in (ksT c /J — H /J) plane get narrower. In addition, increasing 
Hq/J values reduce the dynamical critical temperature while as u increases then the 
F phase region gets wider in the {ksT c / J — Hq/J) plane for all concentration values. 
These observations can also be explained by the well known mechanisms underlying the 
DPT's. In other words, increasing field amplitude facilitates the phase transition due 
to the increasing energy coming from the oscillating magnetic field which tends to align 
the spins in its direction whereas an increasing field frequency gives rise to a growing 
phase delay between the magnetization and magnetic field (i.e. the magnetization 
cannot follow the oscillating magnetic field) and this makes the occurrence of the DPT 
difficult, as a result of this mechanism critical temperature increases. Moreover, the 
low concentration phase diagrams in Fig. (|3J, which appear at small u values such as 
w = 0.1 and 0.5, exhibit a reentrant behavior and the reentrance disappears at a certain 
value of increasing frequency between 0.5 < u < 1.0. For example, as seen from the 
curve with Ho/ J = 0.2, u = 0.5 and c = 0.6, corresponding to the upper right panel Fig. 
01]), as the temperature increases starting from a value in the paramagnetic phase then 
the system undergoes a second order phase transition due to the thermal agitations and 
remains at a ferromagnetic order for a while. As the temperature increases further then 
the system undergoes another second order phase transition from a ferromagnetic to a 
paramagnetic phase. 

An interesting result is that the F + P region disappears in the (k B T c /J — H /J) 
plane for sufficiently weak dilution. In order to focus on this interesting behavior, we 
depict the phase diagrams in the same plane for u = 1.0 with the concentration values 
c = 1.00,0.98,0.96,0.94,0.92 in Fig. (J3J). As seen in Fig. (|3D, as the concentration 
value decreases then F + P region which appears at low temperatures and high field 
amplitude values, gradually shrinks then disappears at a certain concentration value. 
Then we can not see this region in the phase diagrams which are shown in Fig. (jl]) 
except the concentration value c = 1.0. On the other hand, one observes a dynamic 
tricritical point only in a narrow range of site concentration c which disappears at a 
certain concentration value. Namely, as seen in Fig. (j4]) for c = 0.94 we observe a 
reentrant behavior of first order instead of the coexistence region while for c = 0.92 the 
first order reentrance turns into a second order reentrance. 

In order to probe the concentration dependence of the dynamic nature of the system, 
we plot the phase diagrams in (ksT c / J — c) plane for some selected values of Hamiltonian 
parameters in Fig. (J3J). We can see that the results shown in Fig. (j3j) are completely 
consistent with the results presented in Fig. (j3j). As the concentration c of magnetic 
atoms decreases then the ferromagnetic region gets narrower, as expected. In Fig. (j4j), 



we can also observe that there are mainly two reentrance regions which are located for 
Hq/J = 0.2 with u = 0.1, 0.5 and for Hq/J = 1.2, 1.4 with u = 0.5, 1.0, 5.0, respectively. 
The former is nothing but just the reentrant region shown in the upper panels in Fig. (j4]) 
that observed for low concentration values. By taking into account these observations, 
we can conclude that the system may exhibit reentrant behavior for a wide range of 
Hamiltonian parameters as a result of a competition between the interactions that exist 
in the system. Hence, in the paramagnetic phase, the increasing temperature inhibits 
the ability of response of the magnetization to the external field even at low frequency 
region by allowing a transition from paramagnetic phase to ferromagnetic one. 

At a certain c value the system undergoes a phase transition at zero temperature 
which means that the system cannot exhibit an ordered phase below this certain c 
value for the phase diagrams in (fcsT c / J — c) plane. No matter the phase diagrams in 
(ksTc/J — c) plane shows reentrance or not, for concentration values lower than this 
certain concentration value, namely site percolation threshold value (c*), the system 
cannot exhibit an ordered phase at all. In other words, for c > c* the system forms 
an infinite cluster of lattice sites. However, as c gets closer to c* then isolated finite 
clusters appear, since the spin-spin interaction energy become insufficient to compete 
with the magnetic field energy, even in the absence of the thermal energy (i.e. at zero 
temperature). As a result of this, the system cannot exhibit long range ferromagnetic 
order even at zero temperature for c < c*. The variation of c* as a function of the 
magnetic field amplitude for frequencies u = 0.1,0.5,1.0 and 5.0 can be seen in Fig. 
(jl]). As seen in Fig. (Jl]), Hq/J and u parameters have an opposite effect on c*. We 
can say that the general effect of increasing Hq/J is to increase the value of c* whereas 
increasing value of u parameter decreases the value of c*. An important observation 
in Fig. 01]) is the appearance of plateaus for different u values. For u = 0.1, c* value 
is observed between 0.5 < Hq/J < 1.0 which almost does not change its shape with 
varying Hq/J values. As u> value gets higher (e.g u = 5.0 curve on left panel in Fig. 
dl])) then another wider plateau originates within the range 0.0 < Hq/J < 1.0. We can 
see from the right panel in Fig. (jl]) that the same situation holds for a square lattice 
(z — 4). We should state that in the absence of magnetic field {Hq/J = 0), c* values 
are exactly the same values with the static case obtained by the same approximation 
(DA) [33j. 

It is well known that one cannot write an appropriate expression for the free energy 
including the presence of the time dependent oscillating magnetic field. Hence, we 
cannot use the free energy expression to determine the type of the phase transition 
(first or second order). Hence, we simply follow a procedure described briefly below. 
Since the time averaged magnetization over a full cycle of the external magnetic field 
acts as the DOP then we may check the temperature dependence of it. Strictly speaking, 
if the DOP decreases continuously to zero in the vicinity of critical temperature, this 
transition is classified as of the second order whereas if it vanishes discontinuously then 
the transition is assumed to be of the first order. The results of the previously published 
works have pointed out that all transitions from F or P phases to F + P phase (or vice 
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versa) are of the first order transitions [21)1 I2"6] . Hence we shall only investigate the 
transitions between F and P phases, especially for the Hamiltonian parameters which 
exhibit reentrance. The variation of the DOP with the temperature for some selected 
values of the Hamiltonian parameters can be seen in Fig. (j4]). From the left panel in 
Fig. (BJ, we can see that all transitions on the c = 0.6 phase diagram with u = 0.5 
corresponding to the case in the upper right panel in Fig. (J3J) are of the second order. 
The situation is similar for the middle panel of Fig. (J3J) where H /J = 0.2 and u — 0.1 
corresponding to the upper right panel in Fig. (jl]). Finally, we can conclude with the 
transitions shown in the phase diagrams with u) — 1.0 and H /J = 1.4 which are of 
the second order as seen in the right panel of Fig. (Jl]). A detailed investigation for the 
variation of the DOP with the temperature for other Hamiltonian parameters shows 
us that all phase transitions of the site diluted system are of the second order except 
the transitions from F or P phases to F + P phase which occurs at low temperature 
and high field amplitude values and within a very narrow region of the Hamiltonian 
parameter values (e.g. for to — 1.0, c = 0.94 curve as shown in Fig.(T4j)). 

4. Conclusion 

In this work, we have investigated the dynamic nature of the critical phenomena which 
is observed for a site diluted Ising ferromagnet defined on a honeycomb lattice (z = 3) 
driven by an external oscillating magnetic field by means of EFT. We have used a 
Glauber-type stochastic process to describe the time evolution of the system. We have 
given the global phase diagrams, including the reentrant phase transitions. According 
to our calculations, F+P phase disappears for sufficiently weak dilution of lattice sites. 

EFT method takes the standard mean field predictions one step forward by taking 
into account the single site correlations which means that the thermal fluctuations are 
partially considered within the framework of EFT. Although all of the observations 
reported in this work shows that EFT can be successfully applied to such nonequilibrium 
systems in the presence of quenched site disorder, the true nature of the physical 
facts underlying the observations displayed in the present work (especially the origin 
of the coexistence phase) may be further understood with an improved version of the 
present EFT formalism which can be achieved by attempting to consider the multi site 
correlations which originate when expanding the spin identities. We believe that this 
attempt could provide a treatment beyond the present approximation. 

In conclusion, we hope that the results obtained in this work would shed light on 
the further investigations of the dynamic nature of the critical phenomena in disordered 
systems and would be beneficial from both theoretical and experimental points of view. 
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Figure 1. Phase diagrams in (kBT c /J — Ho/ J) plane for a honeycomb lattice with 
some selected values of u and c. The solid (dotted) lines correspond to the second 
(first) order transitions and solid symbols denote the dynamic tricritical points. 
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Figure 2. Phase diagrams for a honeycomb lattice in (ksT c / J — H /J) plane for 
u> = 1.0. The numbers accompanying each curve denote the value of site concentration 
c. The solid (dashed) curves correspond to the second (first) order transitions and 
solid symbols denote the dynamic tricritical points. 
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Figure 3. Phase diagrams for a honeycomb lattice in (kBT c /J — c) plane for some 
selected values of lo and Hq/ J. The numbers accompanying each curve denote the 
value of the amplitude H /J of the oscillating magnetic field. 
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Figure 4. Variation of the site percolation threshold c* of honeycomb and square 
lattices with the amplitude H /J of the oscillating magnetic field. The numbers 
accompanying each curve represent the frequency w of the oscillating field. 
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Figure 5. Variation of the DOP Q for a honeycomb lattice as a function of the 
temperature for some selected values of the Hamiltonian parameters w, Hq/J and c. 



